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Abstract 

We propose a general scheme to construct multiple Lagrangians for com- 
pletely integrable non-linear evolution equations that admit multi-Hamilton- 
ian structure. The recursion operator plays a fundamental role in this con- 
struction. We use a conserved quantity higher /lower than the Hamiltonian 
in the potential part of the new Lagrangian and determine the corresponding 
kinetic terms by generating the appropriate momentum map. 

This leads to some remarkable new developments. We show that non- 
linear evolutionary systems that admit A^-fold first order local Hamiltonian 
structure can be cast into variational form with 2A^— 1 Lagrangians which will 
be local functionals of Clebsch potentials. This number increases to 3N — 2 
when the Miura transformation is invertible. Furthermore we construct a 
new Lagrangian for polytropic gas dynamics in 1 -|- 1 dimensions which is a 
local functional of the physical field variables, namely density and velocity, 
thus dispensing with the necessity of introducing Clebsch potentials entirely. 
This is a consequence of bi-Hamiltonian structure with a compatible pair of 
first and third order Hamiltonian operators derived from Sheftel's recursion 
operator. 

1 Introduction 



In this paper we shall point out a general technique for the construction of 
inequivalent solutions to the inverse problem in the calculus of variations. 
We shall show that completely integrable partial differential equations in 
1-1-1 dimensions that admit multi-Hamiltonian structure can be cast into 
variational form with multiple Lagrangians. It is remarkable that all these 
new Lagrangians can be obtained directly from our present knowledge of 
complete integrability of the evolutionary system without doing any new 
calculations! 



1 



One of the important properties we expect from a completely integrable 
system is multi-Hamiltonian structure. A vector evolutionary system can 
then be cast into Hamiltonian form in more than one way 

{i = 1, 2, n 
K = l,2,..'iV. (1) 
a = —1, 0, 1, oo 

where the variational derivative is denoted by 6k = S/6u'' and J is a matrix 
of differential operators satisfying the properties of a Poisson tensor, namely 
skew-symmetry and Jacobi identity. For integrable systems there exists more 
than one such Hamiltonian operator and Hamiltonian function as the respec- 
tive Hebrew and Greek indices indicate. Then, by the theorem of Magri 
completely integrable systems admit infinitely many conserved Hamiltonian 
functions which are in involution with respect to Poisson brackets defined by 
compatible Hamiltonian operators. 

The essential element in the multi-Hamiltonian approach to integrability 
is the construction of the Hamiltonian operators themselves. Fortunately 
this is a rich subject that can be put to good use. We shall be interested 
in the consequences of multi-Hamiltonian structure on the Lagrangian for- 
mulation of completely integrable evolutionary equations. We shall work in 
the opposite direction to the traditional approach of deriving Hamiltonian 
structure from a Lagrangian. The crucial fact that we shall exploit is the 
relationship between Hamiltonian operators and Dirac brackets [§] for de- 
generate Lagrangian systems which was first pointed out by Macfarlane 
In the case of completely integrable systems we have much more information 
on Hamiltonian structure than Lagrangian and it became clear only recently 
0, 0, 1^ how we can construct multiple Lagrangians for systems that ad- 
mit multi-Hamiltonian structure. We shall now present the general and most 
simple technique for generating these new Lagrangians. 

2 Mult i- Lagrangians 

Evolutionary systems d^) cannot be cast into variational form with a local 
expression for the Lagrangian using the velocity fields alone but require the 
introduction of Clebsch potentials. In 1 + 1-dimensions the general expression 
for Clebsch potentials is given by 

= (2) 
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and in this paper we shall only consider Lagrangians that are local functionals 
of these potentials. In the time-honored way we shall split the Lagrangian 
density for eqs.(|I]) into two 

C = T -V (3) 

that consist of the kinetic and potential pieces respectively. For the first 
Lagrangian density, an enumeration which will become clear presently, the 
kinetic term is always given by 

Ti = 9^k ^\ ^\ (4) 
where gi\^ are constants with det gi\^ 7^ and 

Vi = (5) 

is the Hamiltonian density. We note that the Hamiltonian function that 
appears in (|ip is the space integral of the density. We shall number the con- 
served Hamiltonians by reserving the subscript 1 to the "usual" Hamiltonian 
function but of course there exists conserved quantities such as Casimirs and 
the momentum which are of lower order. In fact, for complete integrability, 
an n-component vector evolutionary system (|lD must admit n infinite series 
of conserved Hamiltonians. We shall denote their densities by 

7Ya;[i] i = l,2,..n; a = -l,...,oo (6) 

and recall that each series starts with a Casimir 

'W-l;^ = 9ikU^ (7) 

which will carry the label minus one. One of these series is distinguished 
in that it contains the "usual" Hamiltonian function which is the one that 
appears in eq.(|ID. For the 2-component systems that we shall discuss in 
this paper these are the Eulerian and Lagrangian series. We note also that 
the two series may coincide up to a relabelling dictated by the recursion 
operator. This is in fact the case for the 7 = 2 case of gas dynamics and 
in most examples of completely integrable dispersive equations except the 
Boussinesq equation. 

The potential part of the Lagrangian does not depend on the velocities 
and from eq.(H) it follows that the Hessian 



det 
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vanishes identically. We have therefore a degenerate Lagrangian system and 
in order to cast it into Hamiltonian form we must use Dirac's theory of 
constraints 0, or the covariant Witten-Zuckerman theory p, ^ of symplectic 
structure. In particular, the first Hamiltonian operator obtained from the 
first Lagrangian is given by 

j{^ = g'^D D = -^ (8) 
ax 

where g'^'^ is the inverse of the coefficients in the kinetic part of the first 
Lagrangian which is non-degenerate. 

The construction of multiple Lagrangians relies on the use of the Lenard 
recursion relation which is implicit in eqs.(|I]) that in the Greek and Hebrew 
indices we have a symmetric matrix 

J'i^6\kiH^] = (9) 

where square brackets denote complete skew-symmetrization and bars enclose 
indices which are excluded in this process. Provided we can invert these 
Hamiltonian operators, we can construct recursion operators 

RX\ = jH':iJtr' (10) 

that map gradients of conserved Hamiltonians into each other (P). 

For the construction of Lagrangians we start with the crucial observation 
that the first Lagrangian is of the form 

A = 0^ - 27^l (11) 

which is manifest from . The original fields that enter into the evolutionary 
system (|I|) are Casimirs which is evident from the subscript minus one. The 
second Lagrangian will be of the same general structure as (|lID if we further 
suppose that eqs. (0) can be written in bi-Hamiltonian form. Thus there will 
exist H2 which is the next conserved Hamiltonian function in the hierarchy 
and the momentum Hq which comes after Casimirs. The higher Lagrangian 
should simply be 

but there is an important refinement that we need to insert here. It is 
not the conserved density but rather the momentum map that enters into 
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the kinetic part of the Lagrangian. The two differ only by total derivatives 
which is irrelevant in the context of conservation laws and therefore generally 
skipped over. However, these divergence terms are of crucial interest as the 
momentum map in the theory of symplectic structure. We shall show that 
given Q*'^ local Hamiltonian structure, the full new Lagrangian is simply given 
by 

Ca = {'Ha-2[{\ + iGa-2[{\)x}(l)l ~ (12) 

where Qa[i] is a functional of the potentials. The coefficient of (j)] above is the 
momentum map and this is the only calculation necessary to find the new 
Lagrangian. 

The fact that it is the momentum map rather than the conserved density 
that plays an important role in the Lagrangian can be seen at the level of 
the first Lagrangian. Now the Casimirs play the role of the momentum map 
and they are used to construct the next higher conserved quantity according 
to the construction of the canonical energy-momentum tensor 

dCi • 1 • 

^0 = "^^T*^^ = = -g' H-i[iiH-i[k] 

which is the momentum. This classical result for Lagrangians linear in the 
velocity can be generalized at each level we have a higher Lagrangian. We 
have 

2Ha-l = g'^[H.a-2[i] + {Qa~2[i])x]'H-l[k] (13) 

ending at the level where a local Lagrangian is no longer possible. In fact the 
validity of this equation is directly related to the existence of the Lagrangian. 
If a check of (|T^) fails for some a, then there exists no local Lagrangian at 

Now we come to an important reservation that our new Lagrangians will 
necessarily carry. The Euler equations that follow from the variation of the 
action with the second Lagrangian will be 

= (14) 

so that the first variation of the second action will certainly be an extremum 
for the original equations of motion (|l|) but the Euler equations (|T4D require 
something weaker, namely linear combinations of functionals in the kernel of 
the recursion operator can be added to the right hand side of the equations 
of motion and the new action will still be an extremum. 
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From this construction it is manifest that for every Hamiltonian function 
in the infinite hierarchy of conserved Hamiltonians that we have for com- 



pletely integrable systems, there exists a degenerate Lagrangian ( [T^ ) that 
yields the equations of motion as its Euler equation up to functionals in 
the kernel of the recursion operator. The number of Lagrangians that can be 
constructed in this way is therefore infinite in number. Given bi-Hamiltonian 
structure we have two local Hamiltonian operators but the Lenard recursion 
operator ( p!OD is non-local. However, the special form of the first Hamiltonian 
operator (H) leads to a local expression for the second Lagrangian in terms 
of Clebsch potentials. But it is clear that the repeated application of the 
recursion operator will require the introduction of non-local terms in higher 
Lagrangians. Strictly speaking, this is not a problem because the original 
Lagrangian is itself non-local in terms of the velocity fields which are 
the original variables. We swept this problem under the rug by introducing 
Clebsch potentials. Higher Lagrangians for evolutionary equations (P can 
be written in local form by introducing potentials for the Clebsch potentials 
themselves! If, however, the equations of motion admit local Hamiltonian 
operators, then our construction guarantees the existence of N Lagrangians 
which are local functionals of the Clebsch potentials. Thus we have 

Theorem 1 A completely integrable system that admits N-fold local first 
order Hamiltonian structure can be given N different variational formula- 
tions with degenerate Lagrangians that are local functionals of the Clebsch 
potentials. 

By a convenient abuse of language we claim that we have a Lagrangian for 
an equation that involves fields when the Lagrangian is in fact only a func- 
tional of the Clebsch potentials for these fields. Then we have the audacity 
to put in by hand the expression for the fields in terms of potentials after 
the variation! This can be at best only a shorthand for the real variational 
principle where we must impose the relationship between the fields and their 
potentials through Lagrange multipliers. 

So far we have been guilty of this abuse ourselves. But now we must say 
that the first Lagrangian is actually 

4f = <j>l, 0^,...) + K{u^ - (15) 

so that upon variation with respect to all the variables (f)^,u^, Aj we get (^, 
Aj = and we arrive at the equations of motion (|1]) expressed in terms of 
the original fields without fudging. 
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Now this obvious observation may seem correct but naive, however, we 
shall now find that it dramatically increases the number of new Lagrangians 
we can construct for integrable systems. 

For every evolutionary equation that admits, say for simplicity, bi-Hamil- 
tonian structure there exists a differential substitution 



u 



(16) 



that brings the second Hamiltonian operator to the canonical form (|^) of 
Darboux. This differential substitution is a Miura transformation. Strictly 
speaking the theorem of Darboux remains unproved in field theory where the 
number of degrees of freedom is infinite but we shall assume it. Miura trans- 
formation works in a direction opposite to the usual action of the recursion 
operator. It leads to Hamiltonian equations 



(17) 



where Hq is the momentum for eqs.(|l|) expressed through (|l^). These are 
modified equations, different from the original equations, but the two sets 
are related by 



(18) 



up to functions in the kernel of some matrix differential operator Oj. 

A comparison of eqs. QI^) and Miura's relation (0) shows us that using 
the differential substitution of Miura we can obtain new Lagrangians for 
nonlinear evolution equations that admit multi-Hamiltonian structure in the 
opposite direction to our earlier construction. Transforming to the variables 
r* and using Clebsch potentials 

r' = (19) 
we can write the classical Lagrangian for the modified system (llT 



u'^=M"'{iP^,il)^^...) 



(20) 



where the labelling of TYq refers to its expression in the original variables 
but these need to be substituted for in terms of r* according to (|T6|) and 
expressed through the potentials (|l^). We note that the Casimirs r* = ipl for 
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the modified system are absent in the polynomial H-aiu) hierarchy. We would 
expect naively that the Euler equations resulting from the first variation of 
the action with the Lagrangian (|20| ) would result in the modified equations 
([T7|). This would indeed be the case if we were to impose the constraint 
between the fields r* and their potentials ip^ as in (|l^) but now using (|19|). 



However, by imposing the constraint through Miura's differential substitution 

4"" = + - (21) 

we obtain a new Lagrangian for the original equations (|I]) in the original 
variables m*. We shall use this construction to derive new Lagrangians, in 
particular for KdV in section |^. It is evident that this construction can be 
extended when there exists multi-Hamiltonian structure but, as we shall find 
in the example of KdV, sometimes it is possible to arrive at local Lagrangians 
using non-local Hamiltonian operators as well. Now we conclude 

Theorem 2 The first Lagrangian of every modified equation obtained 
through a Miura transformation will serve as a new zeroth Lagrangian for 
the original equations of motion provided the constraint between the fields 
and their potentials is imposed through a Miura-type differential substitution. 
For N-fold Hamiltonian structure there exists N — 1 such new Lagrangians. 

Miura transformation is in general not invertible because it is a differential 
substitution. But there exists interesting examples where it reduces to a 
point transformation which is invertible. In that case we can construct — 1 
further Lagrangians. 

We conclude that for an evolutionary system that admits fold first order 
Hamiltonian structure, the number of different variational principles where 
the first variation will be an extremum by virtue of the original equations of 
motion is 2N — 1 and in the case Miura transformation is invertible 3A^ — 2. 
We illustrate this situation for the case of bi-Hamiltonian structure in tables 
^ and 0. The general situation is much more complicated than what these 
tables would lead us to expect. Starting with tri-Hamiltonian structure the 
individual entries in each one of these tables will need to be table by itself 
because there are inequivalent Hamiltonian operators that yield the same 
equations of motion with the same Hamiltonian function. We shall discuss 
this interesting situation in a future publication on the Chaplygin-Born-Infeld 
equation. 
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equations of motion 


Ut = J2SH0 


Ut = J16H1 


J2J1 ^ut = -S^H2 


local Hamiltonian op. 


J2 


Ji 


no 


local Lagrangian 


no 






modified equations 


rt = J2SH0 


n = JiSH, 




local Hamiltonian op. 


J2 = Jl 


Ji 




local Lagrangian 


Co 


no 





Table 1: The hierarchy of local Hamiltonian structures and Lagrangians 
which are local functionals of Clebsch potentials for evolutionary system 
Uf — Ji5Hi — J26HQ where Ji is in the canonical Darboux form. 



equations of motion 




Ut = J2SH0 




local Hamiltonian op. 




J2 




local Lagrangian 




no 




modified equations 


J1J2 = -S^H_i 






local Hamiltonian op. 


no 


J2 = Jl 




local Lagrangian 




Co 





Table 2: When the Miura transformation is invertible we need to include 
an additional column to the left of table 1. 
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3 KdV 



KdV stands as the symbol of completely integrable systems. We think we 
know it, but it turns out to be so rich that there is still new information to 
be learned about it. We recall that KdV 

ut + 6u - u^^^ = (22) 

admits the Kruskal sequence of conserved Hamiltonian densities 

(23) 
(24) 

(25) 

(26) 



^KdV 


= U 






njKdV 
/to 


1 2 

= -u 
2 






'^KdV 


= + 


1 2 

2^^ 




njKdV 
112 


5 4 

= r 


+ 5 M M^. + 


1 2 

2 '^xx 



which are in involution with respect to Poisson brackets defined by two 
Hamiltonian operators 

= D, J2 = -D^ + 2uD + 2Du (27) 

that form a Poisson pencil. By introducing the potential 

M = 0, (28) 

KdV can be cast into variational form with two Lagrangians 

^Kdv ^ H5r0,_2<^^ (29) 
jO^d^ = (^H!^^v^^^^^^^_^^Kdv (30) 

which consist of the classical Lagrangian and the second Lagrangian re- 
spectively.Q Here we observe that both ( [29| ) and ( pO[ ) are examples of our 
general expression ([T2|) for higher Lagrangians. 

The second application of Lenard's recursion operator to Ji results in a 
third Hamiltonian operator which is non-local so we cannot continue to gen- 
erate higher Lagrangians. But we can use Theorem 2 to generate new lower 

^Notc that there is an error in potential term of the second Lagrangian in |^. The 
results that follow are stated correctly. 
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Lagrangians for KdV. For this purpose we note that in both Lagrangians 



p9| ) and (pOD we should have added the constraint X{u — (jy^) and written 
the full Lagrangian. But following the convenient abuse of language we did 
not do so because it was manifest. It is, however, necessary to write the full 
Lagrangian in the case of lower Lagrangians. 

According to our general construction of lower Lagrangians we first recall 
the original Miura transformation 

u = r"^ + (31) 

that brings J2 to the canonical form of Ji in the variable r. The equation of 
motion for r is mKdV which is different from (^) but under the substitution 
([3l| ) we have Miura's result 

Ut + 6m Ua, - Ua:xx = {D + 2r) (vf + 6r^r^ - r^^^'j = (32) 

so that, on shell, if mKdV is satisfied then so is KdV. Now we can introduce 
the Clebsch potential for the modified field variable 

r = i'x (33) 
and write the first Lagrangian for mKdV 

CT'"''' = ^xi^t + K''' (34) 

in a straight-forward manner. But now enforcing the constraint in the full 
Lagrangian through the Miura transformation 

^ W full ^ ^^K,V ^ _ ^2 _ (35) 

we shall arrive at a new Lagrangian for KdV because the Euler equation that 
comes from the first variation of this action will be satisfied by virtue of (^2]). 
Unlike ([3^) which is a higher Lagrangian, (|35D is a lower Lagrangian in the 
sense of the action of the recursion operator on the equations of motion in 
the resulting Euler equation. 

And the saga of KdV continues! We consider the third Hamiltonian 
operator for KdV 

J3 = R^Ji (36) 

which is nonlocal but the relationship between differential substitutions and 
Hamiltonian structures of KdV |^| enables us to construct another new 
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local Lagrangian for KdV. For this purpose we recall that the differential 

substitution ^ 

r = aq-\ h (37) 

q 2q 

which transforms mKdV into twice modified KdV 

= (fc - ^ + ^ - 2aV^ (38) 
is a Miura transformation for (|36|). This can best be seen by the expression 



J3 = ^{q'^D + Dq^) - q^D~\^ (39) 

for the third non-local Hamiltonian operator for KdV in terms of twice mod- 
ified variable q. We recall that J3 is fifth order in u. We have the Miura 
relation 

2 _ i ^ 'l^ ^ ^ 



rt + 6r - r^^^ = [a - — + -D] 

qt- [qxx-^ + 2a2g3 



2q q 



(40) 







between modified and twice modified KdV's. Introducing the potential for 
the twice modified variable g = we have 

2 2 

u = $(Xx, X.., X...) = ^ - ^ + 2ax.. + a\l + 2ae + ^ (41) 

a: Az Ax 

in terms of the original field u. The first Lagrangian for twice modified KdV 
is simply 

Cr''''' = X.Xt + H^r , (42) 

and therefore the second lower Lagrangian for KdV is given by 

full ^ ^rn.KdV ^ ^[^^ _ ^(^^^ ^^^^ (43) 

which provides another illustration of (|2T|). This process can be continued. 

We note that an alternative to the Clebsch potential for KdV is the 
Schwartzian which was pointed out by Schiff [jl0| . We shall postpone consid- 
eration of Schwartzian potentials to future work. 
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4 Polytropic gas dynamics 



The simplest examples for applying our construction of multi-Lagrangians 
consist of quasi-linear second order hyperbolic equations that Dubrovin and 
Novikov ||lT| have called equations of hydrodynamic type. The distinguished 
example in this set consists of the Eulerian equations of polytropic gas dy- 
namics in 1 + 1 dimensions 



pt + up^ + pu^ = 
ut + uu^ + p^^'^p^ = 



(44) 



and in particular for 7 = —1 we have the case of Chaplygin gas, or Born- 
Infeld equation that was recently shown to have a string theory antecedent 
]r2| . This system can be cast into quadri-Hamiltonian form [^. For the 
Chaplygin-Born-Infeld case the complete Hamiltonian structure can be found 
in and its symmetries were given in In the following we shall use 
the labelling = p and = u. 

First we have three local Hamiltonian structures of first order 



16 



■h 



pD 




Du + {-f 



Dp 

-2)uD 



(7 



-2)Du + uD 
-^D + Dp^-^ 



u p D + D u p 



D 



1(7 



v 



+ 



D 

i(7 



-2W 



2" ^ 7 



2)n^ + :^P' 

,7-1 



,7- 



D 



D 



u p' 



D 



D u 



(45) 



(46) 



(47) 

which form a Poisson pencil J = Ji + C1J2 + C2J3 with Ci,C2 constants, 
i.e. these Hamiltonian operators are compatible. In eq.(^) is the Pauli 
matrix and this is the canonical Darboux form of first order Hamiltonian 
operators. The equations of polytropic gas dynamics admit two infinite hi- 
erarchies of conserved Hamiltonians which are in involution with respect to 
Poisson brackets defined by all three of these Hamiltonian operators. In the 
first set, which is called Eulerian [TH], the Hamiltonian densities are given by 



p 



(48) 
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= up (49) 
«f = r""+T^f (50) 

= r'^f+TH^f' (51) 



where (^) is the momentum, (|50D is the famihar Hamiltonian function, the 
Casimir is in P8| ) and the rest consist of higher Hamiltonians. Therefore, 
the Euler series is the distinguished one in the terminology of section 0. The 
second series 

1-L^_^ = u (53) 
< = l{^-2y + -^p^-' (54) 

= - 2y + p-y-' (55) 

D 7 — 1 



is the Lagrangian series which starts with the Casimir (|5B|). Note that for 
7 = 2 this series is no longer polynomial as logarithms will enter and the 
same remark holds for integer and half-integer values of 7 in both series. 

Finally, we note that the recursion operator R2 = J2Ji^ can be used to 
write infinitely many Hamiltonian operators by letting it to act n times on 
J I. However, in general none of these operators will be local. In particular 
we note that 

= J3 {Ji)-' 7^ {Rif, J3 ^ J2Jl^J2 (57) 

except in the case of shallow water waves where 7 = 2 which admits extension 
to integrable dispersive equations. 

Next, there is a third order Hamiltonian operator [117 1 which was obtained 
from Sheftel's remarkable recursion operator 

= DU-^ a^D (58) 
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where 



U = { _A_V2 ^ ) (59) 



7-2 



p' U 



which is only compatible with Jq. Higher conserved Hamiltonians start with 
the density [|l|], pf 

U^-P^ ^Px 

which is part of the Eulerian series. There is also a Lagrangian series starting 
with 

m''^ = —T^^, (61) 

U^-P^ ^Px 

and both form new infinite hierarchies of conservation laws. 

We will be interested in the Lagrangian formulation of the equations 
of polytropic gas dynamics (^^ that correspond to all these Hamiltonian 
structures. Introducing the Clebsch potentials 

M = V^x, P = V^x (62) 
we have the first Lagrangian representation for this system 

C[ = n'l.^t + n^m - m^{^,, ^.) (63) 

but using the recursion operators J2 Jf ^ and J3 Jf ^ we find two further 
Lagrangians 

CI = n^^t + ni^t-2ni{^,,^,) (64) 
CI = n^iJt + n^vt-mfiv^.ij,) (65) 

which are local functionals of the Clebsch potentials. The Lagrangian ob- 
tained through the action of the recursion operator J4 Jf^ is the most inter- 
esting one. Because J4 is a third order operator, the fourth Lagrangian 

d = n'X^''^u, + n'X^'^^pt-2n^,{^.,iJ,) (66) 

^7 pxUt - U^pt ^ ff.^. 

^4 = — zzn - 2 p 67 

- P^ ^pI 

is local in the velocity fields. This is a general property of bi-Hamiltonian 
structure with a pair of first and third order Hamiltonian operators. Here 
we find a remarkable situation in that the number of Lagrangians that we 
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can construct by repeated application of Sheftel's recursion operator J^Ji^ 
is infinite in number. All of these Lagrangians will be local in the original 
field variables p and u. 

Now we come to lower Lagrangians that will arise from Miura transfor- 
mations. The Miura transformations that bring the Hamiltonian operators 
(56) and (^7]) to the Darboux form of (|^ ) are point transformations for 
equations of hydrodynamic type. Dubrovin and Novikov had pointed out 
that first order Hamiltonian operators for equations of hydrodynamic type 
are given by 

r' = g''D-g^-Tl^u: (68) 

where Qik are the components of a Riemannian metric which is fiat by virtue 
of the Jacobi identities. The Miura transformation provides manifestly flat 
coordinates for this metric. For example from ( ^6] ) we find the fiat metric 

dsl = _ _ 1)2^2 - {-f -l)udpdu + p du'^\ (69) 

and it can be verified that the Miura transformation 

p = rp u = (r^^^ + p'^'^M (70) 

7 — 1 ^ ^ 

brings it into the manifestly fiat form 2dr dp. In these variables we find the 
first modified equations of gas dynamics 

n + (r^"' + P^-y. + irp^-^p, = (71) 

7-1 

Pt + 'jpr'^~'^r^ H rir'^'^ + P^^^)Px = 

7 - 1 

and linear combinations of these equations with variable coefficients give 
eqs.(0) of gas dynamics. Introducing the potentials 

r = Xx, P = v^ (72) 

we have the Lagrangian 

£2^"" = Xxvt + v.xt - + \[u- ^^"^^f"' ) + ^ (P - XxVx) (73) 

where is the momentum (^) expressed in terms of the potentials x 

V. Transforming to the first modified variables r, p we get Ji, J2 = Ji and J3 
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defining the tri-Hamiltonian structure of eqs.(|7T]). Now tliere is a new lower 
Lagrangian tliat we can construct from the recursion operator JiJ2~^. We 
find 



Ji 



[r^-i + (7 - 2)p^-i] AD 
V +L'[(7-2)r^-i+p^-i] A 



1-7) [pr^-2^L) + Dpr^-'^/\] 



A 



1 



(74) 



(7- l)(r^-i -pT-i)2 

where the labelling of the variables is in the order r and p. The new La- 
grangian is given by 



7/u// _ Xx'^t — fJxXt 



CI'; 



Xi - vi 



7-1 , 



+ 0- (P - XxVx) (75) 



where the momenta do not belong to the polynomial series of conserved 
Hamiltonians. However, we can identify the lower momenta from this La- 
grangian 



3-7 



where ± refers to Eulerian and Lagrangian series as well as the roots of the 
quadratic equation. 

We now turn to the third Hamiltonian structure (P7| ) defined by the fiat 
metric 



dsl 



K7-i)^ 



[(7 - 1)2^2 - Ap-y- 

^ ^(7 - Ifu^ + 4p^"^] dp du + up du^] (76) 



2(7-1) 
2dq dw 



and the coordinate transformation that brings it to the manifestly fiat form 
is given by 



(7 - lyu^ - Ap 



o7-l 



7-3 
2(1-7) 



(77) 
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w = I — ^e^de (78) 




1 (7 - l)u + 2/)(^-i)/^ I 

2 (7 - 1)m - 2p(7-i)/2 j 

where, in general, the last integral cannot be done in closed form. For some 
specific values of 7 the integral (^) is elementary as in the notable case 
of Chaplygin-Born-Infeld. But this paper is devoted to the general case of 
polytropic gas dynamics and we shall not consider inverting (|77D, ([78|) to 
obtain u,p as functions of q and w. We shall only remark that after this 
inversion we can obtain two more new Lagrangians. 

The Lagrangians (pHD, ( |51D and (^) for polytropic gas dynamics are ex- 



amples illustrating the general expression ([12|) for higher Lagrangians. For 
equations of hydrodynamic type there is no dispersion and hence Q van- 
ishes identically. We have given only two (ffSD, ([75|) of the four lower La- 
grangians because the integral (^) must be carried out before we arrive at 
the second modified equations of gas dynamics which will lead to two fur- 
ther new Lagrangians. Certainly the Lagrangian ( |67|) which is derived from 
bi-Hamiltonian structure with a first and third order operators according to 
([I^ ) is the most remarkable one because this is the first time it has been 
possible to write a Lagrangian for polytropic gas dynamics that is local in 
the original field variables, namely the density and velocity. Furthermore it 
is only the first element in an infinite series of such Lagrangians. 



5 Kaup-Boussinesq system 

Gas dynamics with 7 = 2 governs the behavior of long waves in shallow 
water. From the point of view of complete integrability it is a remarkable 
case, because in this case we find several completely integrable dispersive 
generalizations of eqs.(^^. Most prominent among them is the well-known 
Kaup-Boussinesq system E^ 



Ut 




pt = (up + e'^u^^^^ (79) 
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which admits tri-Hamiltonian structure. The first Hamiltonian structure is 
given by the Hamiltonian operator (^Sj) and 



jKBq 



D 



\Du 



\uD UpD + Dp)+e^D'' 



(80) 



where denotes the principal value integral, is the second Hamiltonian 
operator for the Kaup-Boussinesq system. In the limit e — > this Hamilto- 
nian operator reduces to p6|) with 7 = 2. The recursion operator is given 
by 

\u + \u.,D-^ 1 



TjlKBq 
K2 



^ , (81) 

and there is a third local Hamiltonian operator obtained by the action of the 
recursion operator J^'^'^ = {RI ^^'^YJq as in the 7 = 2 case of gas dynamics. 
The conserved Hamiltonians in the Eulerian and Lagrangian series are 



KBq 



HI 



^KBq 



P 

u p 

1 

2 
1 

2 

1/1 Son 1q 49 

-U'p + + -p3 + ^4^2^ 

-^'^{\pul + ^uu^p^ + pI + ^u\l) 



pu^ + p^ + e^u u^^ 



pu^ + 3p^M — e^{Auxpx + 3mu^.) 



(82) 
(83) 

(84) 
(85) 
(86) 



and the degeneracy in the 7 = 2 case of gas dynamics is repeated in its 
dispersive generalization. In particular, the Lagrangian and Eulerian series 
coincide apart from a relabelling 



^KBq{E) 



KBq{E) 



U 



^KBq{L) 
-j^KBqiL) 



(87) 



njKBq{E) 



njKBq{L) 



that is dictated by the recursion operator. 
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With the aid of the Clebsch potentials 



ti = V^x, P = i^x (88) 

we obtain 

^KBg ^ T^^^^Yt + H^i'^^t - 2ni^%(p^, Vx, Vxx, i^xx, •••) (89) 

for the first Lagrangian. Using the technique we have presented in section 
we shall now construct higher Lagrangians. These three local Hamiltonian 
structures enable us to construct two new Lagrangians 

C^^l = [Hq^'^ + e^Lpxxx)'^t + H^l'^i^t - 2H2^'^{<^x, i^x, ^xx, 'ipxx, •••) (90) 



and 



n^'"' + ^Vxo:.) i't - m^''\v,, ...) (91) 



for the Kaup-Boussinesq system. The determination of Qp-ii] is according to 
eq.(|l2D with /3 = 2, 3 and [2] = [1] — 1 because of the relabelling difference (|87D 



between the Lagrangian and Eulerian series. Note that the momentum map 
which is the coefficient of (pt in (|90|) is exactly the same as the momentum 
in front of ipt in (PT|). The reason for this goes back to the degeneration 
of the Eulerian and Lagrangian series into one and the fact that it is the 
momentum map that is the important element in the general construction 
([T^). In the dispersionless limit the Lagrangians (^), (0), (piP reduce to 



the gas dynamics Lagrangians (p3), (|6^) and (p5f) with 7 = 2 



6 Kaup-Broer System 

There is another completely integrable dispersive version of the 7 = 2 case 
of gas dynamics which is the Kaup-Broer system The triangular 

invertible differential substitution 

p = r] + €u^ (92) 

transforms the Kaup-Boussinesq system into the Kaup-Broer system 

ut = uux + r]r, + euxx 

Vt = {m)x - ^Vxx (93) 
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which also has three local Hamiltonian structures For the Kaup-Broer 
system the conserved Hamiltonians in the Eulerian series are given by 

K""' = ur] (94) 
HP'- = \[u'v + v'-'2evu,] (95) 

n^^'' = ^[u^T] + SuT]'^ + 6er]uu^ - ie'^u^T]^], (96) 
n^'' = \u'v + luW + lv' + eilv'u.-u'v.) (97) 

which can be obtained from (^)-(^) through the substitution (P^). The 
first Hamiltonian operator for the Kaup-Broer system is given by ( ^5]) and 
the second Hamiltonian operator 

Jr'=(i rf n2 l^n^'n\] (98) 
1 yi'uD-eD^ l{r]D + Dr]) J ^ ^ 



can be obtained from (|80|) of the Kaup-Boussinesq system using the substi- 
tution (^). 

For Kaup-Broer system we introduce the potentials 

7] = Wx, 1p = W + Eifix (99) 

and arrive at the first Lagrangian 

= + n'^2'-wt - 2nf'''{w,, w^x, ^,x, ...) (100) 

but now we can derive two further Lagrangians using the recursion opera- 
tor obtained from the Hamiltonian operators (^) and (P5|). Following our 
procedure of section Q we find the second Lagrangian 

C^^r ^ ^^KBr _ 2ew,,)vt + U'^f'^Wt - 27^f ^^(u;,, w,,, ip^x, ...) (101) 

which is the same as the Lagrangian of Kaup-Boussinesq system ( PTf ) subject 
to the differential substitution (0). Similarly we find 

= «^'^ + 2£2«;_)y,, (102) 

as the third Lagrangian for the Kaup-Broer equations (0). As in the case 
of Kaup-Boussinesq, these Lagrangians reduce to 7 = 2 gas dynamics La- 
grangians in the dispersionless limit. In the Kaup-Broer Lagrangians we find 
another example of the general formula (|12D for Lagrangians. 
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7 Nonlinear Shrodinger equation 



We shall consider the nonlinear Shrodinger equation in the 2-component real 
version 

,2 



V 



V 

Y 

Vt = {v^)x, (103) 

which is a reaction-diffusion system. Again this reduces to the 7 = 2 case of 
gas dynamics in the dispersionless limit. This version of NLS can be obtained 
by another triangular differential substitution 

u = v + erix/ri (104) 

from the Kaup-Broer system. 

NLS has the same first local Hamiltonian structure (^) as in the case of 
Kaup-Boussinesq or Kaup-Broer systems. Once again the second Hamilto- 
nian operator for NLS can be found by the transformation ( |104| ) from the 
second Hamiltonian operator (|98|) of the Kaup-Broer system. Thus for the 
2-component real version of NLS the second Hamiltonian operator is given 
by 

l[iv-%.D + D{7]- 



J. 



NLS 



D 



2 



and the conserved Hamiltonians are 



njNLS 


= V 




= V 


njNLS 


= VT] 


'^NLS 


-\ 


njNLS 

I-L2 


1 

2 



rjv + 7] 



_2^r 



njNLS 
^3 



riv + 3vri + e 



3 9 9 X o 1 

4 4 




(105) 

(106) 
(107) 
(108) 

(109) 
(110) 
(111) 
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which forms an infinite sequence combining both Eulerian and Lagrangian 
series according to (p^). 

In order to construct the Lagrangians for NLS we introduce the potentials 



V 



z = If — elnWx 

and the first Lagrangian 

r^LS _ njNLS ^ ,njNLS^„ On-lNLS/^,, ^ 



;ii2) 



;ii3) 



is the classical result. Once again we shall use the techniques of section | 
to construct higher Lagrangians with the recursion operator obtained from 
). We obtain two higher Lagrangians for NLS 



([T05|) and 



rNLS nj 

L,r, = /T| 



NLS , 
+ 



Wt 



u;i-2H^^^K,z.,...) (114) 



and 



rNLS 



■ Wt 

3 (2^x5 Ulx, •••)• 



;ii5) 



2^NLS 



that are local functionals of the potentials. Here again, in the dispersionless 
limit we find the 7 = 2 gas dynamics Lagrangians. The remarkable strength 
of the general expression (|1^) for new Lagrangians is manifest. 



8 Boussinesq Equation 

In order to discuss the bi-Hamiltonian structure and the Lagrangians for the 
Boussinesq equation in a unified framework we first turn to its dispersionless 
limit. For polytropic gas dynamics we had 

pt = ipu)x, = I y + ) ^^^^"^ 



with its first nontrivial commuting fiow 

^ .7-2 



Py = Ux, Uy = (117) 
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both of which reduce to a second order quasi-hnear wave equation [|T^ . If we 
express Boussinesq equation in the form 



or 



Pyy 



Py 



1 1 

P Pxx 



p 







e'pxx 



118) 
;il9) 



as a first order evolutionary system and compare its dispersionless hmit to 
polytropic gas dynamics, we find that it corresponds to the commuting flow 
for 7 = 4. The completely integrable dispersive equation 



Pt 

Ut 



pu — 2e 



u 



2 1 



2 +3^ 



(120) 



X 



is the commuting flow to the Boussinesq equation. 

This system admits bi-Hamiltonian structure p5[ with the Hamiltonian 
operators (EH) and 



J 



f pD + Dp- Se'^D^ 

3Du — 2ux 
\ 



3uD + 2ux 



\ 



8{p^D + Dp^) + 8e^D^ 
-e^5ipD^ + D''p)-3{pxxD + Dpxx)] J 



;i21) 

which are compatible. The conserved Hamiltonian densities for the Boussi- 
nesq system are given by 

(122) 
(123) 







njE 


= pu, 


nf 


1 ■ 


4 . 


njE 

ri2 


1 

28 
+2 



4„2 

XX 



2pu' + -p^ + e\Qppl + Aul) + 4£V; 
14 7 

—pu^ + -p% + Ue'^{2uul + 4pVxMx + ^uppl) (124) 

(125) 



in the Eulerian sequence and we have also 



(126) 
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< = u^ + \p' + e'pl (127) 
= + -p^M - £'u(4pp,, + 3p^) + — eVxPxx, (128) 

^2 = \u^ + \p''u' + ^p'' + e\^p'pl + Au'pl + ?,2(mp,u, + Sp'ul) 
4,136 2 2 248 4 128 16 , 96 ^ x / n 

+e (— p p^^ - — p^ + — MMxxPxx + yM^^Pxx + yPM^^)(129) 

g 2 592 3 64 2 ^ 64 8 2 

-\-e \Olpp^^^ Y^P^^ ~^ ~^^xxx) "I" Pxxxx 

in the Lagrangian sequence. The Hamiltonian function of Boussinesq system 
with the first order Hamihonian operator in Darboux form ( ^5]) is \H^. 
We note that the system (|119|) for the Boussinesq equation differs from all 
dispersive integrable examples we encountered earlier in that its familiar 
Hamiltonian function ( |128| ) is in the Lagrangian sequence. This is because 
Boussinesq equation is the family of commuting flows to the regular gas 
dynamics hierarchy. The first commuting higher flow for the Boussinesq 
system ( |120| ) has the Hamiltonian function ( |124| ) in the Eulerian series. 
By introducing potentials 

M = V?a:, P = i^x (130) 

we can obtain two local Lagrangian densities for the Boussinesq system. First 
we have the classical Lagrangian 

for Boussinesq system and its first nontrivial commuting flow ( |12(]| ). The 
second Lagrangians are given by 

C2^^^ = {Hq ^^'^ - 4£: Vxxx) V^y + [^0 - {i^xi^xx)x + 4^ " 'H{ 

(133) 

£f ^) = 7=4_4^2^^^^)^^^ [^L 7=4_5^2(^^^^^)^^4^4^^^^^J^^_2^i^ 7=4 

(134) 

according to the general construction of Lagrangians in ([121). Here we see 
also that the Lagrangian for the commuting flow is obtained by flipping 
the Hamiltonian functions between the Lagrangian and Eulerian series while 
keeping the momenta fixed. In section we had constructed Lagrangians 
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for gas dynamics using the Hamiltonians from the Eulerian series in the 
potential part of the Lagrangian. The general formula (^2|) can readily be 
used to construct Lagrangians for the commuting flow ( |119| ) by this simple 
flip in the potential. 



9 Conclusion 

This is the first time it has been possible to write a Lagrangian for polytropic 
gas dynamics that is local in the original field variables, namely the density 



and velocity. It is a result of the general expression (12) that serves to identify 
immediately multi-Lagrangians for completely integrable systems. What is 
even more remarkable is that this is only the first element in an infinite series 
of such local Lagrangians for polytropic gas dynamics. 

It is worth emphasizing again that the scheme we have presented in sec- 
tion^ is a universal one for the construction of multi-Lagrangians appropriate 
to evolutionary systems. The expressions (|T2D and ( pT] ) for Lagrangians of 
completely integrable systems has general validity. We note that {\L2^ with 
a = 1 is true even in the case of non-integrable equations, provided the equa- 
tions are presented in the form of conservation laws and the system admits 
one further conserved quantity, namely the Hamiltonian. We have discussed 
in detail the higher Lagrangians for the completely integrable non-linear evo- 
lution equations of polytropic gas dynamics, Kaup-Boussinesq, Kaup-Broer, 
NLS and Boussinesq equations all of which bear out the universal applicabil- 
ity of ([T2D in the construction of higher Lagrangians. We have also presented 
the lower Lagrangians (pi]) fully for KdV and partially for gas dynamics ow- 
ing to the difficulty of writing the second modified variables in closed form. 

The invariance group of these multi-Lagrangians and their Noether cur- 
rents should prove to be of interest in discovering new hidden symmetries 
of fluid mechanics. We did not discuss this important issue here. Recently 



Jackiw and co-authors |^ have used hidden symmetries in the classical La- 
grangian for fluid mechanics to construct very interesting field theory models 
of fluid mechanics. Multi-Lagrangians may prove to be of interest in this 
connection also. 
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